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Abstract 

We have revealed evident errors in the statistical analysis, performed by Klapdor- 
Kleingrothaus et al in a recently published paper [J to establish a limit on the stability 
of electron with respect to the decay into ^ + 7. The performed reestimation of the sen- 
sitivity of the experimental setups to the 256 keV gamma emitted in the hypothetical 
electron decay, has shown that the limits on the electron stability and charge noncon- 
servation parameter e^-y presented in [l], have been overestimated by at least a factor 
of 5. 

PACS: 13.35.-r; 13.40.Hq; 14.60.Cd; ll.30.-j 
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Incongruity of the analysis is evident already from the comparison of the sensitivity estimate 
with the "lcr" method and the results of the maximum likelihood and x 2 - analysis, the latter 
two give the result by a factor of 3-5 better. The most evident manifestation of errors in 
the analysis can bee seen in Table 7 of Ref.|l] (called hereafter the KKKT article): 

1. The best fit of the ANG2 data set contains 89.444±63.058 events^] in the peak corre- 
sponding to ~256 keV (m e /2) 7 from hypothetical electron decay e — > v + 7 for the 
ME casc|. As follows from the text of the article the error in this number is cited 
at 68% c.l. (it is called an indication of a signal on a 1.4 u c.L). The corresponding 
upper limit on the number of events in the peak, A, already at 50% c.L, should be 
higher than the central value of 89 in most practical cases of any almost-symmetrical 
X 2 - profile shapes. Instead of the above it is claimed to be only 38 events at 68% c.l. 

1 We are reproducing the number of significant digits, following the original text of the KKKT article. 
2 See the original text of the KKKT article pQ to explain the abbreviations, we are also citing a number 
of values from the KKKT article without going too much into details. 
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2. The upper limit on the number of events in the 256 keV peak obtained for ANG4 set, 
is 4.789 events at 68% c.l. This value is lower than la statistical error on the content 
of a single bin in Fig. 8d, while the FWHM peak width is about 20 bins and the mean 
bin content is > 50 events. 

These evident errors, together with the fact that an "indication of a signal on 1.4 a" for 
1.94 xlO 26 yr is excluded both by the results of the Borexino [2] and DAMA [3 J collaborations 
with more than 90% probability, convinced us to have a closer look at the results of KKKT. 
The results from the KKKT article are summarized in Table Q] (the data are taken from 
Tables 5 and 7 of the KKKT article; for the sake of simplicity we are citing only the ME case 
at 68% c.l.). These are the results of the analysis of experimental data with the la method 
("ler" column) and with the standard least square procedure ("x 2 " column). One can see 
incompatibility of the values in N pea k and A subcolumns, the most evident discrepancies are 
described above. 

In order to check the achievable sensitivity when looking for the gaussian shape on the 
linear background, we have applied a toy model consisting of a gaussian peak superimposed 
on the linear background. For the toy model we used the background level for the corre- 
sponding data set from the KKKT article, the region of analysis was set to 100 keV with the 
bin width of 0.36 keV, and the la width of the gaussian peak corresponded to the Doppler- 
broadened line of 3.25 keV. The number of the events that can be eliminated at la level 
was defined for a large number of samples by using the MC method. For each sample a set 
of randomly distributed events was simulated with the fixed mean number of events in the 
gaussian peak, then we fitted it with linear+Gauss analytical shape. The sensitivity at 68% 
confidence level in this approach corresponds to the mean number of events in the peak for 
which the x 2 value increases by A% 2 = 1. The results are presented in "MC" column of 
Tabled 

Comparing the A values in MC column to the KKKT la estimate one can see that 
the toy model gives the same level of sensitivity, confirming the KKKT estimation with 
the la method. Nevertheless, the values of A obtained by KKKT with the x 2 method are 
significantly (by 2-10 times) lower than their own estimates with the la method. 

There is no description of the % 2 - profile analysis in the text of the KKKT article. It is 
also not clear how many free parameters were used, what is the precise n.d.f, and whether the 
result depends on the lower and upper limits of the analysis region. The number of the events 
in the peak ("peak area") seem to have been taken directly from the minimization program 
together with the error on this number. In principle, the correct limits with acceptable 
precision can be reproduced by using these data. The results are presented in Tabled The 
data of A are obtained assuming the normal shape for the % 2 - profile with the central value 
and variation taken from the third column of Table 7 of the KKKT article (these values 
are reproduced in N pea k column of Table [J) . Our limits have been calculated by using the 
Bayesian approach (see i.e. [4], the prior knowledge in our case is the restriction of the 
positively defined effect). The limits on the life-time recalculated by using these values are 
presented in the next column. Only the ME case for 68% c.l. is shown, analysis for 90% c.l. 
and the AE cases can be performed in the similar way with the same conclusions. 

As it is seen from our estimations the best limit that can be obtained (r = 3.9 x 10 25 yr at 
68% c.l. for ANG4 setup) is very close to that already existed for the HPGe (r = 3.7 x 10 25 
yr at 68% c.l. [5]). The most surprising fact is that it practically coincides with the estimate 
of the sensitivity obtained by KKKT themselves with the la method. As a result, the life- 
time limits in the KKKT article from the x 2 analysis are by a factor of 2 till 10 stronger 
than the estimated sensitivity, and by a factor 3 till 5 higher than the values that can be 
obtained from the values presented as the "Peak area"H 

3 The same is true for the limits on the charge- nonconservation parameter £^„ 7 , derived from the upper 
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T 
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("peak area") 


events 


years 


ANG1 


49 


5.8 x 10 24 


-38. 187 ±51.077 


13.216 


2.146 x 10 2b 


ANG2 


61 


3.3 x 10 2b 


89.444 ± 63.058 


38.354 


5.285 x 10 2b 


ANG3 


64 


2.2 x 10 25 


-38.301 ± 67.374 


13.216 


10.76 x 10 25 


ANG4 


46 


2.0 x 10 25 


-76.249 ±47.401 


4.789 


19.33 x 10 25 


ANG5 


73 


3.0 x 10 25 


-33.273 ±75.947 


14.343 


15.69 x 10 25 



Table 1: Data from the KKKT article (68% c.L; ME case only). A is an upper limit on the 
number of events in the peak from the hypothetical electron decay, r is the corresponding 
life-time limit. 
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Table 2: Reestimation of the life-time limits (68% c.L; ME case only). 



In the proper analysis one should take the signal from all electrons (the AE case) into 
account, not only from the outer shells, that can obviously change the sensitivity. The limit 
for the AE case is inferior to the existing for the HPGe, so there is no need in a more detailed 
analysis. 

Some words should be written on the model used to describe the underlying background. 
The quality of the fit for 4 sets is bad (x 2 — 390/280, this value of % 2 from the formal 
statistical point of view, rejects the model with a very high probability) and only for the 
ANG2 set it has an acceptable quality (\ 2 — 280/280). If the data are obtained under the 
same conditions (which seems to apply at least to 4 detectors of Setup 1), then the model 
should give a statistically compatible description for all sets. The quantitative comparison of 
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the data sets can be performed using Fischer's F-distribution ^ = F{p, v, v) as a significance 
test, where v is a number of the degrees of freedom and p is a confidence level (see i.e.|6j). 
Solving equation F(p, 280, 280) = 390/280 with respect to p, one obtains the statistical 
probability of the data set with lower x 2 '- P — 0.003. This definitely points on the systematics 
problem with the data set ANG50. 

The linearity of the background has not been justified in the KKKT article, moreover, 
all the compatible sets contain the statistically evident hole in the background just in the 
place where the effect is searched for. 

limits on the electron life-time. Using the best established limit on the electron life-time (T e ^„ 7 > 4.6 x 10 26 
y, 68% c.l. and the upper limit on the photon mass (m 7 < 7 X 10 -19 eV [!]), from formula e 2 _^„ 7 = 

(j^) 5 T 6xl ° a 3 one can obtain the best restriction on the parameter e 2 ^„ 7 < 0.23 X 10 -98 . 

4 or, if we assume that the linear model of the underlying background is valid, the situation is inverted, 
and in this case the data for all detectors except ANG2 should have systematic problems. The question 
what the real situation is should be addressed to the KKKT authors. 
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Conclusions 



Our analysis of the data presented by KKKT in [T] has shown that the upper limit on the 
electron decay is overestimated by at least a factor of 5. The statistical analysis of the 
KKKT contains evident errors. Moreover, one of the presented data sets (containing an 
"indication of a signal" on 1.4 a) is statistically inconsistent with 4 others, pointing out on 
possible systematic problems with the experimental data. The model used in [1] to fit the 
underlying background has the same problems of statistical incompatibility. 

The best limit that can be obtained by using the KKKT data is comparable to the one 
established previously for the HPGe. The KKKT restriction for charge nonconservation 
£ e~>u~t < 0-86 x 10~ 98 is not valid either, since it is based on the overestimated electron 
life-time. Instead of the above the restriction < 2.3 x 10~" at 90% c.l. can be 

calculated from the modern best limit on the electron life-time r e ^ vl > 4.6 x 10 26 yr (90% 
c.l.) established by the Borexino collaboration. 
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